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I Introduction 



The main goal of this paper is to perform spherical averages over the hyperbolic coordinates in 
Minkowski spacetime by means of direct integration. We work with Lorentz invariant (or "radial" ) 
functions, defined on M^'", which depend only on the distance to the light-cone — {x^^Y ^ i^f'- 
(.r" is the time; coordinate and x stands for the spatial coordinates). In particular, we will compute 
their Fourier transforms; 

F{k) := / /(r?^.a;''a;^) exp(-27rifc^a;'') , (1) 

where the metric 'q^j^v is of the form 

r?^,=diag(l, -!,...,-!) (2) 

(We set the 27r in the exponential to avoid factors of 2-k appearing in front of the integral. Here it 
is just for convenience but this strategy becomes crucial when integrating over infinite-dimensional 
spaces.^) 

This type of integral is ubiquitous in Quantum Field Theory (QFT), since the fields can be 
decomposed in terms of their Fourier components, and the correlation functions depend only on 
the radial distance s^. Although QFT is defined in Minkowski space, a common procedure is to 
perform a Wick rotation at an early stage (this procedure can be found in standard textbooks, 
e.g., Ref. 2) in which all integrals over spacetime are reduced to Euclidean integrals. 

There are several obvious disadvantages to relying on the Wick rotation procedure, in which 
the time coordinate t is rotated to imaginary values t — > —it, and Minkowski spacetime is trans- 
formed into Euclidean space. First, the special structure of Minkowski spacetime in which there 
is a preferred direction (i.e., time direction) is lost. Second, Wick rotation might not be a valid 
procedure in all circumstances, e.g., in cases where the arc at infinity does not vanish. 

The motivation for this work is to provide tools for performing calculations in QFT while 
working directly in Minkowski space. 

In analogy with Euclidean space, we introduce a (pseudo-) spherical coordinate atlas, i.e., one of 
the angles is hyperbolic. Integrating over all the angular variables leaves us with a one-dimensional 
radial integral. The advantage of this procedure is that any possible light-cone singularities of the 
integrand are mapped to point singularities in the radial variable and therefore easier to deal with. 

In this paper we will not specify the function space for which the integrals are convergent. We 
will assume that only such functions have been chosen. 

Below, we will discuss the cases n = 1 and n = 2 in detail. For the case of general n we only 
quote the results. 

After preparing an earlier version of this paper we found an article by Codelupi'^ which seems to 
be largely unknown. Codelupi studied Fourier transforms of Lorentz invariant functions /(s), using 
related methods, and found the same results. To find the case of arbitrary n he proved a recursion 
relation, that relates Fourier transform of n -|- 2 spatial dimensions to the Fourier transform of n 
spatial dimensions. In Section IV we adapt Codelupi's elegant method to our case. 

II Case Ri'i 

II. 1 Pseudospherical coordinate atlas for M^'^ 

A characteristic feature of manifolds with an indefinite metric is that a global spherical coordinate 
system does not exist.'* To cover R^'^ we use four patches and four different parametrizations. 
The parametrizations coincide on the boundary between the domains, i.e., on the light-cone. (See 
Fig. 1). 

Denote by t,x the global Cartesian coordinate system with distance = — x"^. We 
parametrize the different patches as follows: 
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Figure 1: Standard division of R^-^. I is the forward light-cone, III the backward hght-cone, and 
II and IV the spacelike domains. 



• Patch I: 

t = s cosh X = s sinh V', 

where s G [0, +oo[ and t/j g] — oo, +oo[. 

Volume element: dt A dx = s ds A dip. 
Lino clement: s"^ = — x^. 

• Patch II: 

t = is sinh ip, x = is cosh tjj, 

where s e]ioo, iO] and e] — oo, +oo[. 
Volume clement: dt A dx = s ds A dtp. 
Line clement: = — x"^. 

• Patch III: 

t = s cosh tp, X = s sinh tp, 

where s G [0, — oo[ and tp g] — oo, +oo[. 
Volume element: dt A dx = sds Adip. 
Line elem( 

• Patch IV: 

t = is sinhf/), a; = is cosh ■^j 



where s g] — ico, iO] and xp g] — oo, -|-oo[. 
Volume element: dt A dx = s ds A dip. 



The limits of integration in each patch are chosen to yield a positive result when the volume element 
is integrated over a small, finite volume. 

II. 2 Fourier transform of radial functions 

We integrate the function /(s^) separately in the timelike and spacelike domain for the cases of 
timelikc and spacelikc momenta. To simplify the calculation, we set kx = when the momentum is 
timelike, and kt = when the momentum is spacelike. This can always be achieved with a Lorentz 
transformation and is not a restriction on the results. 

• Patch I+III: 
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(i) Timelike momentum: = 0, 

Ii+iii{kt) = / dt dx f (t^ — x"^) exp{—2'Kiktt) 

Ji+iii 

f-OO f + OO 

= / ds s /(s^) / dip exp{—2'KiktScoshtp) 

Jo J-oo 

r — oo p+oo 

+ / ds s f[s^) / dijj exp{—2'iriktScoship) 

Jo J-oo 
poo P+OO 

= 2 dso So /{sq) / dip 008(2^ kt s cosh ip) 

Jo J-oo 

= -2-K dso So f{sl) No{2TTktSo), (3) 
Jo 

where Nq is a Bessel function of zeroth order and 



so = Vt^-x^- (4) 

In the last line we have used formula 3.868(2) from Ref. 5, after a change of variable to 
x = e'l'. 

(ii) Spacelike momentum: fct = 0, 

Ii+iii{kx) = / dtdxf(t^ — x^) exjp{—2Trikxx) 
Ji+iii 

POO P+OO 

= / rfs s /(s^) / dtp exp{—2'KikxS sinh. ip) 

Jo J-oo 

p — oo r+oo 

+ / ds s f[s^) / dtp exjp{—2mkxssinhip) 

Jo J-oo 

POO p-\-oc 

= / dso So /(sq) / dtp cos(27rA;xSo sinhf/;) 

Jo J-oo 

POO 

= 4/ dsoso f{sl) KQ{27rkxSo), (5) 
Jo 

where Kq is a Bessel function of zeroth order and sq as before. In the last line we have 
used formula 3.868(4) from Ref. 5 , after a change of variable to a; = e'^. 

Patch II+IV: 

(i) Timelike momentum: kx = 0, 

Iii+ivikt) = / dtdxf(t^ — x^) exp{—2'Kiktt) 

J II+IV 

iO p+oo 

ds s f[s^) / dip eyip{2TrktSsinhip) 

loo J-oo 

iO p+oo 



/lu p-\-oo 
ds s /(s^) / dtp exp(27rfctssinh^) 
-ioo J-oo 
POO p+oo 

= I dsi Si f[s^) I dip cos{2'KktSism]iip) 

Jo J-oo 

POO 

= 4 / dsi Si f{sl) Ko{2TrktSi), (6) 
Jo 



where Kq is a Bessel function of zeroth order and 

si = - t2. (7) 

In the last Hne we have used formula 3.868(4) from Ref. 5, after a change of variable 

X = e^. 

(ii) Spacelike momentum: fct = 0, 

Iii+iv{kx) = / dtdxf(t^ — x^) exp{—27Tikxx) 
Jii+iv 

= / ds s f[s^) / d-tp eyip{2Tr kxS cosh ijj) 

Jico J —oo 

/id p+oo 
ds s f[s^) / exp(27rfcxScosh^) 
-ioo J —oo 

poo f+OO 

= 2 dsi Si /(si) / dip cos(2'KkxSicos]i'tp) 

Jo J -oo 

POO 

= -27r / dsisi f{sX) No{2T^kxSx), (8) 
Jo 



where A'^o is a Bessel function of zeroth order and Si as before. In the last line we have 
used formula 3.868(2) from Ref. 5, after a change of variable to a; = . 



In summary 



poo POO 

I{kt) = -2w dsoso f{sl) No{2nktSQ)+A dsi si f {si) Kq {2-KktSi) 
Jo Jo 

POO POO 

I{kx) = 4 / dsoSof{sl) Kq {2TTkxSo) - 2n / dsiSif{sl) No{2TTkxSi) 
Jo Jo 



II. 3 Example: Fourier transform of a Gaussian 

In this section we apply the results from above to a specific test function and show that the correct 
answer is obtained. Whenever necessary, we define the integral dx f (x) as lim^^ci Jq dxe^*^^ f{x) 
With this proviso, we can directly compute the Fourier transform of a "Gaussian" (where the quo- 
tation marks remind us that the exponential is imaginary): 



dtdxe'^' )e-2vrifctt ^ ^g-^TT fc,_ 
We can now check that our Fourier transform integrals give the same result, 

POO 

7i(fcj) = 4/ drre-^^+''"'^K^{2'Krkt) 
Jo 



1 1 



:r(i + |)r(i-f)e'T^i^ 1 ^ 



POO 

I2 {kt) = -27r / drre-^'-'^''^N^{2Trrk 
Jo 



kt ^fl^i sin(^) 



e 



5.U2(e-i); ^■''ni + v) 5.|V2(e-i) 
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where we have used 6.631(2,3) from Rcf. 5. Each integral has a pole as — » 0; the W function 
with negative first argument has a simple pole in i^, and there is an inverse sine in v in the second 
integral. However, as wc will see, this singularity exactly cancels between the two integrals. Thus 
we may, in fact, take v Q in the sum. First we convert the M function to W functions, 

2'2 \2'' 2'2 \ ~'~2-' 2'2 

using formula 9.233(1) in Ref. 5. Then we make use of the following identities (9.234(1,2) and 
9.235 in Ref. 5): 

Wi .{z) = i+j,{z) + ^W i^iz), 

2' 2 Z ^> 2 Z ^1 2 

W i.{z) = i+^{z), 

~2'2 ^ "'2 

"■2 

Expanding the F functions and 1/ sinz = 1/z + z/6 + ©(z^), we can explicitly verify that the 
singularity cancels: 

2» 9? 

/i+/2 = — + - + 0{\) = + 0{\) 

V V 

We have verified numerically that the constant term is, in fact, 

• 2/2 

T I T — ITT K 

li + 12 = ire 

III Case 2 

III.l Pseudospherical coordinate atlas for R^'^ 

Global coordinate system: t, x, y, distance =t^ — — y^. 

• Patch I: 

f = scosh^, a; = s sinh i/) cos ^, ?/ = ssinh'^sin^, 

where s G [0, +oo[, ijj £] — oo, +oo[ and 6 E [— 7r/2, tt/2]. 

To avoid problems when sinh'!/; switches sign at 0, the integral over ijj needs to be split up 
into two integrals: ijj € [0, +oo] and ip £ [0, — oo]. 
Volume clement: dt A dx A dy = sinh^ Adi/j A dO. 
Line element: =t^ — — y^. 

• Patch II: 

t = issinht/), a; = is coshT/! cos6', y = iscosh'i/'sin^, 

where s e [iO,ioo[, V — oo, +oo[ and 9 e [—■ 7r/2,7r/2]. 
Volume element: dt Adx Ady = is^ cosh^ A dip A dO. 
Line element: = — x^ — y^ . 

• Patch III: 

t = scosh'^, a; = s sinh i/) cos ^, ?/ = s sinh sin 0, 

where s G [0, — oo[, ^ — co, +oo[ and 9 G [— 7''/2, 7r/2]. 

To avoid problems when sinh?/> switches sign at 0, the integral over ^ needs to be split up 
into two integrals: g] + oo, 0] and '0 g] — oo, 0]. 
Volume element: dt Adx Ady = sinh^ Adij) A dO. 
Line element: =t^ — x^ — y^ . 
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• Patch IV: 



t = is sinhtlj, x = is coship cos 6, y = iscoshipsinO, 

where s G] — ioo, iO], -0 g] — oo, +oo[ and 9 £ [— 7r/2, n/2]. 
Volume element: dt A dx A dy = is^ coshtp ds A dip A d9. 
Line element: s^ = i'' — x^ — y^ . 

The limits of integration in each patch are chosen to yield a positive result when the volume 
element is integrated over a small, finite volume. 



III. 2 Fourier transform of radial functions 

As before, wc integrate the function /(-s^) separately in the timelike and spacelike domain for the 
cases of timelike and spacelike momenta. To simplify the calculation, we set = ky = when 
the momentum is timelike, and kt = ky = when the momentum is spacelike. This can always be 
achieved with a Lorentz transformation and is not a restriction on the results. 

• Patch I+III: 

(i) Timelike momentum: kx = ky = 0, 



U+III 



(kt) 



i+iii 







dtdxdy fif' — x^ — y^) exp{—2mktt) 

" POO 

ds s^ f[s'^) / rf-i/) sinh ■^exp(—27rifcts cosh 



+ / dip sinhipexp{—2'!riktSC0ship) 
Jo 

r—oo 

+7r / ds s^ f{s^ 
Jo 







+ 



dip smhipexjp{—2'iriktSC0ship) 

r 

/ dip smhipexp{— 2wikts cosh ip) 

J — oo 



ds S^ /(i 





2 

h Jo 



cxp(— 2 TT ;/.:/, s) 
ikts 



ds s'f{s') 



ikts 



dso sq /(sq) sin(27rfctSo) 



where 



So = \Jt'^ — x'^ — y^. 
The angular integrals have been computed as follows: 



(9) 
(10) 



/•oo — oo 

/ dil) sinh ^exp(—io cosh V') + / dip sinh exp(—m cosh -0) 
Jo Jo 

r+oo—ie 

m.2 I dip sinh7/;exp(— iacosh^) 

2 



lim: 



lim 

la 



exp(— iacosh(oo — ie)) + exp(— iacosh(0 — ie)) 



= — exp(— ia) , 
la 

where a = 2Trkts, and in the last line cosh(V' ± ie) = cosh^ ± iesinh^ has been used. 
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(ii) Spacelike momentum: kf = ky = 0, 

Ii+iii{kx) = \ dt dx dy f {t'^ - x'^ - y'^) exp{-2mkxx) 

= / ds f(^s^) / d6 dip sm]iipexjp{—2mkxSsinhtpcos9) 

Jo J -IT 12 L ^0 



+ 



+ 



/ dxp sinh'!/'cxp(— 27riA:2:Ssinh'0cos6 
Jo 

p—oo r^/2 r pO 

/ ds f[s^) / dO rfV' sinh^exp(— 27rifcxSsinh^cos^) 

Jo J-v/2 iJoo 



rO 

+ dip sinhip exp{—2TrikxS sinhip cos6) 

J — oo 

Jo '^xS Jo 



exp(— 27rfcxs) 



2 f°° 

TT / '^■^0 «o /(sq) exp(-27rfca;So) . 

"^2: Jo 



(11) 



The angular integrals have been computed as follows (a = 2'Kkxs): 

f7r/2 



/ ' 

J--K/2 



POO p — OC 

/ dip sinhip exp(— msinh'i/'cos^) + / dip sinhip exp(— iasinh^/; cos( 
Jo Jo 

poo f'^/2 

4 / dip sinh^ / d9 cos(asinh^cos^) 
JO Jo 

POO 

■-2iT I dip sinh')/; Jo (asinh^/i) 
Jo 



277 



exp(— a). 



To get to the third line we have used formula 3.753(2) from Ref. 5, after a change of variable 
to a; = cos^?. The last line is obtained using 6.554(1) from Ref. 5, after a change of variable 
to y = sinh'^. 

Patch II+IV: 

(i) Timelike momentum: kx = ky = 0, 



I 



II+IV 



(kt) 



/ dtdxdy fif' - - y"^) cxp{-2TTiktt) 

J II+IV 

pioo r'^/2 p+oo 

i / /(s^) / d9 dip cosh^ exp(27rA;tSsinh^) 

JiO J-tt/2 J -oo 



/lU f^/^ r+oo 

ds f(s^) / dO dip coship exTp{2'KktS sinhip) 

-ioo J—'ir/2 J — oo 

ds' s' f{—s' J / coshf/; exp(— 27riA:fs'sinhf/;) 



—IT / ds' s' /(—a' ) / dip cosh^ exp(— 27riA;ts' sinh^) 

J oo J — oo 



0. 



(12) 
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The angular integrals have been computed as follows (with a = 2'Kkts'): 

f + CX) 



dtp cosh^exp(— iasinh^) 

-oo 



/+00— ie 
dip coshV'exp(— zasinh^) 

+00— ze 

exp(— iasinh^) 



= lim — 

e^O a 



= lim — 

e^O a 

= 0. 



exp(— iooo) exp(— eaoo) — exp(iaoo) exp(— eaoo) 



(ii) Spacelike momentum: kt = 0, 

hi+iv = / dtdxdy f (t^ — — y^) exp{—2Trikxx) 
Jii+iv 

pioo 



= I 



( / ds f[s^) / d6 dtp cosh^ exp(27rA;xS cosh V' cos ^) 

JiO J-Tr/2 J -00 

/iO r+oci /'7r/2 

dss^fi^s^^ I dtp cosh tp / d0 exp(27rA;xS cosh ■i/' cos 0) 

-ioo J — CO J—7r/2 

POO roo 

= 4 / dsi si /(— si^) / d9 dtp cosh^ exp(27rifcxSi cosh ^ cos ^) 
Jo Jo Jo 

/•oo foo 

+4 / dsi f[si) / d6 dtp cosh^ exp(—27riA;j:Si cosh ^ cos ^) 
Jo Jo Jo 

pQC r^/2 roc 

= 8 dsi sf f{sf) / d0 dtp cosh')/; cos(27rA;a;Si) 
Jo Jo Jo 



/•oo 

■ / dsisif{sl) COs(27rfca;Si) , 
Jo 



(13) 



where 



Si = \J x^ -\-y^ — t"^- 
The angular integral has been computed as follows (with a = 2'KkxS\): 

f7r/2 



roo r~ I ^ 

I dtp cosh^ / d9 cos(acosh^ cos^) 
Jo Jo 

= — / dtp cosh tp Jo (a cosh ^) 
2 io 



= — cos(a), 
a 

where in the second line formula 3.715(19) from Ref. 5, and in the last line formula 
6.554(3) from Ref. 5 have been used, after a change of variable to a; = cosh^'. 



In summary 
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2 f°° 

I{kt) = -— / dso8Qf{sl) sin(27rfctSo) 
i^t Jo 

2 2 
I{kx) = T dso So f{sl) eyip{-2-jrkxSo) + — ds^ si f (sf) cos{2'!rka:Si) . 

f^'.r Jo If I- Jo 



IV Case Ri'" 

In this section we follow Codelupi's derivation^ of the 1 + n dimensional case. The idea is to derive 

a recursion relation between the Fourier transform in n and n + 2 spatial dimensions, and then use 
the explicit expressions found before to construct the general case. 

IV. 1 Recursion relation 

Define the radius in n spatial dimensions as follows: 



i=l 



Formally, the function /(s) = /{Vf^ — r^) always looks the same, independent of the number of 
spatial dimensions. Suppose we have spaces with spatial dimensions n = 1 to n = m. For each of 
these spaces exists a transform 

/•OO 

fco) = / dr xn{r, k) G(r, ko), (14) 
Jo 

where (proof given in the appendix) 

Xn{r, k) = "^T^ j^^/^-i -^71/2-1 (27rrfc) (15) 

and ^ 

G{r,ko)= I dtf(y/t^-rA exp{-2TTikot) . (16) 

J — OO 

But formally, G{r, ko) looks the same for all cases, for example, 

G{r,ko)= dkXm{k,r)F^"'\k,ko), (17) 
Jo 

assuming that in the inverse Fourier transform, the angular contribution can also be integrated 
out. 

Substitute in the line before, 

POO POO 

F^"\k,ko) = / drxn{r,k) / du Xm{u,r) F^"'\u,ko) 
Jo Jo 

POO POO 

= / duF(^\u,ko) dr Xn{r,k) Xm{u,r). (18) 
Jo Jo 

We can explicitly evaluate the second integral using formula 6.575(1) from Red. 5. The result is 

drxn{r,k)xm{u,r) = ^^u{v?-k'^)^ ^Q{u-k) (19) 
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where h = (m — n)/2, T{x) Euler's Gamma function and Q the step function. Equation (|18|) now 
takes on the form 

F(-Hk,ko) = f^J^ duF(-+^'^Hu,ko)u {u^-k')''-\ (20) 

Considering the special case h ^ 1, i.e. m = n + 2, and taking the derivative with respect to k of 
both sides of this equation leads to the recursion formula 

i^("+2)(fc,fco) - --^T |:^^"^(fc,^o). (21) 

2nk ok 

So the problem is solved, at least in principle, once we find the explicit expressions for n = 1 and 
n = 2. But as we will show in the next section (again following Codelupi^), the recursion relation 
above will also allow us to find explicit formulas for the case of general n. 

IV. 2 Explicit expressions for M^ " 

Let us define 



la = Jkl-B, 



The recursion relation Eq.(^l|) can be rewritten in terms of Zq and ^i. 

To find expressions for general n consider the following two cases 
(a) n-even. 



s, 



{-l)^2n / dso /(so) 4^ (2^50^0) ■ (24) 



^0^ 



The last equation is proved be iteratively applying the derivative to F^"^^ (Iq) (result of Section 
III. 2) and using formula 8.472(2) from Ref. 5. Similarly, we find 



n+l 



F(")(Zi) - 4/ dsa fiso) ^ K^{2^soh) 

n 

-2ti / dsi f{si) ^^ (27rsi/i) , (25) 



7 2 ^ 

where the formulae 8.472(2) and 8.486(13) from Ref. 5 have been used. 
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(b) n-odd, 



/ 1 



V27r/o 



and 



= (-l)'^27r / dso /(so) -i^ (27rso/o) 
Jo l^^ 

n + l 

n-1 Z'"" S ^ 

+(-1) — 4 / dsi /(si) 4^ Jfn^ (27rsi/o) 
Jo / 2 ^ 

i-co -2- 

F(")(Zi) = 4/ dso f{so) ^ K^{27TSoh) 
Jo / 2 ^ 



-277 / dsi /(si) ^ (27rsi/i) ., 

Jo ^ 



(26) 



(27) 



using the results of Section II. 2. and, again, formulas 8.472(2) and 8.486(13) from Ref. 5. 

For even n, the Fourier transform with a timelike momentum has no contribution from the space- 
like region of spacetime. 

We can summarize both cases in the following formulas, now valid for arbitrary n: 



-2tt / dso /(so) -Iht 
+Jr.-i (27rsoZo) sini tt 



n — 1 



A^n-i (27rsoZo) cos^TT ^ 



f(")(Zi) 



n- 1^ ^ 

/■°° s~ / n-1 

+4 / dsi /(si) 4rr -fi^i^ (27rsi;o)cos tt— — 
^0 ' ^ 

/■°° s~ 
4 / ciso /(so) 4rr (27rso;i) 
^0 ^1 ' ' 

n + l 

-27r / dsi /(si) 4rr (27rsi«i) . 

^0 / 2 ^ 
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Appendix : Derivation of Eg . (|T51) 



Equation (jT^) can be derived in several different ways. For direct integration see e.g., Ref. 6, 
Chap. 4. We follow here a derivation presented in Ref. 7. 
We define the Fourier transform on R" by 

/+00 r+oc 
... dxi . . .dxn f{xi, . . . ,X2) exp{-27Ti {kiXi + . . . + knXn)) , (Al) 

-oo J ~oo 

(A2) 



or short 



F{Pj = j (Tx f{r) exp(-27rifc. 
and the inverse Fourier transform by 

f{r) = j d"k F(tj cxp(2TTik-r^ . 

One property we will need is 

A2/(f) = ^ = -47r2 / ci"fc k^ F(k] expf 2nik ■ 



(A3) 



(A4) 



Now assume we have a radial function, and we work in spherical coordinates. We would like to 
write the Fourier transformation as 



F(fc)= / dr f{r)xn{r,k) 
Jo 



and the inverse Fourier transform as 

fir) 



dk F(k) Xn{k,r)dk, 



(A5) 



(A6) 



where the x's contain the integration over the compact n ~ 1 angular coordinates. From Eq.( A4) 
we also know 

/•OC 

A2/(r) = -47r2 / dk k^ Xn{k, r) F{k). 
Jo 



(A7) 



Calculating A /(r) (now starting with Eq.(A6)), 



d'fjr) 
dxl. 



dk Fik) ^ Xn{k,r), 
oxL 



where 



yields 



\r / or r'^ or'^ 



A' fir) = 



d'fjr) 



72 — 1 Q 

dk F{k) ( — — — xn{k, r) + -g^ Xn 



But with Eq.dATD 



dk F{k) 



n — 1 d 

TT^ Xn{k, r) H — Xn{k, r) + 'ii:'^ k^ Xn{k , r) 



0. 



(A8) 



(A9) 



This equation is valid for arbitrary F{k), so the expression in brackets has to be zero. This ODE 
(in r) has the general solution (see Ref. 8, p. 146) 



Xn(fc,r) = A„(fc)ri-"/2Zp(27rfcr) 



(AlO) 
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where A^^k) is determined by the initial conditions, Z is a Bessel function of order p, and p = 
± (1 — n/2). Computing the inverse Fourier transform expUcitly in the cases of n = 1 and n = 2, 
determines p to be n/2 — 1. 
To find An, consider /(r) at r = 0, 

/•OO 

/(0)= / dk F{k) xn{k,0). 
Jo 



From Eq.( AlO ) we have 

Xn{k,0) = limA„(fc)ri-"/V„/2_i(27rrfc) 



r^O 



An{k) 



Anik) 



r(n/2) 
(7rfc)"/2-i 



T{n/2) 



where 9.1.7 from Ref. 9 has been used. 

But according to the definition of the inverse Fourier transform 



/(O) = / rffc"F(fc) 



r(l + n/2) 



/ dk F{k). 
Jo 



The factor in front of the integral is the volume of the unit n — 1-sphere. Equating both expressions 
for /(O), which are valid for arbitrary F{k), yields 

A„(fc) = 2^fc"/2, 

and therefore as the final result 

Xn{k, r) = 2^fc"/2ri-"/2 J„/2_i (2^rfc) . (All) 
Because of the symmetry of the transformation, 

X„(r, k) = 27rr"/2fci-"/2 J„/2_i (2^rfc) . (A12) 
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